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@ Measurement of thermal conductivity and specific heat. 

@ A method and apparatus for determining both the thermal 
conductivity, k, and specific heat, c p , of a fluid of interest are 
disclosed. An embodiment uses proximately positioned resis- 
tive heater 126 and thermal sensor 122, 124 coupled by the fluid 
of interest. A pulse of electrical energy is applied to the heater 
of a level and duration such that both a transient change and a 
substantially steady-state temperature occur in the sensor. The 
k of the fluid of interest is determined based upon a known 
relation between the sensor output and k at steady-state 
sensor temperature; and c p of the fluid of interest is determined 
based on a known relation among k. the rate of change of the 
sensor output during a transient temperature change in the 
sensor and c P . 
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MEASUREMENT OF THERMAL CONDUCTIVITY AND SPECIFIC HEAT 



CROSS REFERENCE TO RELATED APPLICATIONS 5 

Reference is made to two related applications 
published under numbers and (agents ref- 
erences B2012698 and B2012738, corresponding to 
U.S. Patent applications Serial Nos 07/21101 14 and 10 
07/21 1200) filed of even date and in the name of the 
common applicant. 



BACKGROUND OF THE INVENTION 15 



1 . Field of the Invention 

The present invention relates to the measurement 
of certain physical properties of fluids and, more 20 
particularly, to the determination of both the specific 
heat and thermal conductivity of gases. In a 
preferred embodiment a trapped gas sample trans- 
mits steady state and transient responses to input 
energy of limited duration which can be measured 25 
electrically as by extracting the influence of the input 
energy in the form of measurable change in 
temperature of an appropriate sensor in contact with 
the gas of interest. 

30 

2. Prior Art 

In the prior art the traditional approach to 
determining specific heat, c p , has been via calorimetry 
using reversible step increases of energy fed to a 
thermally isolated or adiabatic system. Such devices 35 
are bulky, slow and cumbersome. Little progress has 
been made toward the automation of a rapid method 
to make this determination. 

With respect to measuring thermal conductivity in 
fluids various types of detectors have been used. 40 
This includes resistance bridge type sensors. One 
such device is described in U.S. Patent 4,735,082 in 
which thermal conductivity is detected using a 
Wheatstone bridge technique in which a filament in 
one diagonal of the bridge is placed or positioned in 45 
a cavity through which the sample gas of interest is 
passed. The filament is used to introduce a series of 
amounts of thermal energy into the fluid of interest at 
alternating levels by varying the input voltage which, 
are, in turn, detected at the other diagonal as voltage 50 
difference signals. Integration of the changes of the 
value of the successive stream of signals yields a 
signal indicative of the heat dissipation through the 
fluid, and thus, the thermal conductivity of the fluid. 

Further to the measurement of thermally induced 55 
changes in electrical resistance, as will be discussed 
in greater detail below, especially with reference to 
prior art Figures 1-5, recently very small and very 
accurate "microbridge* semiconductor chip sensors 
have been described in which etched semiconduc- 60 
tor "microbridges" are used as condition or flow 
sensors. Such sensors might include, for example, a 
pair of thin film sensors around a thin film heater. 



Semiconductor chip sensors of the class described 
are treated in a more detailed manner in one or more 
of patents such as 4,478,076, 4,478,077, 4,501,144, 
4,651,564 and 4,683,159, all of common assignee 
with the present invention. 

It is apparent, however, that it has been necessary 
to address the measurement of specific heat c P , and 
thermal conductance, k, of a fluid of interest with 
separate and distinct devices. Not only is this quite 
expensive, it also has other drawbacks. For 
example, the necessity of separate instruments to 
determine specific heat and thermal conductivity 
may not allow the data consistency and accuracy 
needed for useful fluid process stream (gas or liquid) 
characterization because the required degree of 
correlation may not be present. 



SUMMARY OF THE INVENTION 

The present invention overcomes many disadvant- 
ages associated with the determination of both 
specific heat, c p , and thermal conductivity, k, by 
providing simple techniques which allow accurate 
determination of both properties in a sample of 
interest using a single sensing system. The present 
invention contemplates generating an energy or 
temperature pulse in one or more heater elements 
disposed in and closely coupled to the fluid medium 
(gas or liquid) of interest. Characteristic values of k 
and c P of the fluid of interest then cause correspond- 
ing changes in the time variable temperature 
response of the heater to the pulse. Under relatively 
static sample flow conditions this, in turn, induces 
corresponding changes in the time-variable re- 
sponse of one or more temperature responsive 
sensor coupled to the heater principally via the fluid 
medium of interest. 

The thermal pulse of a source need be only of 
sufficient duration that the heater achieves a 
substantially steady-state temperature for a short 
time. This pulse produces both steady-state and 
transient conditions at the sensor. Thermal conduc- 
tivity, k, and specific heat, c Pl can be sensed within 
the same sensed thermal pulse by using the 
steady-state temperature plateau to determine k 
which is then used with the rate of change of 
temperature in the transient condition to determine 
c P . 



BRIEF DESCRIPTION OF THE DRAWINGS 



Figures 1, 2, and 3 are different views of a 
prior art embodiment of a microbridge flow 
sensor. 

Figures 4 and 5 are typical circuits for use 
with the sensors of Figures 1-3. 

Figure 6 is a schematic representation of 
sensor time/temperature response curves ac- 
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cording to a heater pulse. 

Figures 7a, 7b, and 7c, represent several 
heater/sensor configurations of microbridge 
systems in accordance with the invention. 

Figure 8 is a scanning-electron-microscope 
(SEM) photo of the microstructure of a typical 
microbridge sensor. 

Figure 9 is a partial schematic and block 
diagram of a circuit for use with a sensor as 
depicted in Figure 7(b) in accordance with the 
invention. 

Figure 9a is a more detailed circuit schematic 
with reference to Figure 7c. 

Figure 10 is a schematic block diagram of the 
system of the invention including calibration 
and use functions. 

Figure 11 is a scope trace representing the 
temperature signal rise versus time, for the 
configuration of Figure 7(c) in response to a 
heater pulse for dry air at atmospheric press- 
ure, 

Figure 12 is a graphical representation of the 
temperature signal rise versus time, for the 
configuration of Figure 7{c) in response to the 
heater pulse for various gases at atmospheric 
pressure as indicated. 

Figure 13 is a graphical representation of 
thermal conductivity determination based on 
the bridge output of Figure 9(a). 

Figure 14 is a theoretical graphical represen- 
tation of sensor heat-up time versus pressure 
for several gases using the sensor configura- 
tion of Figure 7b. 

Figure 15 is similar to Figure 14 based on data 
taken by a sensor of the type depicted in Figure 
7(b) calculated in accordance with the inven- 
tion. 

Figure 16 is a graphical representation of 
sensor heat-up time versus pressure for several 
gases using the sensor configuration of Figure 
7c. 

Figure 17 is a graphical representation of 
sensor cooling time versus pressure for several 
gases using the sensor configuration of Figure 
7c. 



Detailed Description 

The present invention, then, is directed to a 
system which enables both the determination of 
specific heat. c p , and thermal conductivity, k. The 
system utilizes a thermal pulse approach which is 
based on generating an energy or temperature pulse 
in a heater, which is coupled to a sensor primarily by 
the fluid medium (gas or liquid) of interest. Both 
quantities can be determined from a single pulse. 

Thermal conductivity and specific heat of each 
fluid of interest produce characteristic transient and 
steady-state temperature reactions in a proximate 
sensor as exemplified in Figure 6. 

In the preferred implementation, specific tempera- 
tures, as Ti and T2 in Figure 6, are selected as 
'marker" points with respect to the sensor. These 
marker points are used to reference the determina- 



tion of the time periods, as ti-t2, required to 
achieve the corresponding temperature rise(s) or 
fall(s) in the sensor(s) between the marker points. 
As will be discussed, the sensor or sensors are 
5 located in predetermined spaced relation to the 
heater or heaters, but preferably physically separ- 
ated therefrom so that the proximate influence of the 
solid heater material(s) is reduced and the coupling 
of the heater with the sensor or sensors by the fluid 

10 of interest is relatively enhanced. 

The preferred embodiments of the approach of 
the invention contemplate disposing spaced micro- 
spec sized heating and sensing elements in a 
relatively static (zero flow) sample of the fluid of 

15 interest. The microsensor system or "microbridge" 
system, as it will be referred to herein, though not 
limiting, is presently preferred for several reasons. 
The system is extremely fast reacting, is very 
accurate, very sensitive because of its advantageous 

20 coupling to the fluid of interest and small and 
adaptable to a variety of configurations. 

The microbridge semiconductor chip sensor con- 
templated, for example, in certain embodiments 
preferred for the invention may resemble the form of 

25 one or more of the microbridge systems illustrated 
in the patents identified above. Such a system is 
exemplified by Figures 1-5 taken from Patent 
4,501,144. A discussion of that example will now be 
presented as it will be helpful in understanding the 

30 present invention. While the present discussion is 
believed sufficient, to the extent necessary, any 
additional material contained in the microbridge 
related patents cited is deemed to be incorporated 
herein by reference. 

35 The illustrated embodiment of Figures 1-5 con- 
templates a pair of thin film temperature sensors 22 
and 24, a thin film heater 26 and a base 20 supporting 
the sensors and heater out of contact with the base. 
Sensors 22 and 24 are disposed on opposite sides of 

40 heater 26. Body 20 is a semiconductor, preferably 
silicon, chosen because of its adaptability to preci- 
sion etching techniques and ease of electronic chip 
producibility. The embodiment includes two identical 
temperature sensing resistor grids 22 and 24 acting 

45 as the thin film heat sensors and a centrally located 
heater resistor grid 26 acting as the thin film heater. 

Sensors 22 and 24 and heater 26 may be 
fabricated of any suitable, stable metal or alloy film. 
In Figure 8, the metal used was a nickel-iron alloy 

50 sometimes referred to as permalloy, with a composi- 
tion of 80 percent nickel and 20 percent iron. The 
sensor and heater grids are encapsulated in a thin 
film of dielectric, typically comprising layers 28 and 
29 and preferably silicon nitride, S13N4, to form thin 

55 film members. In'the embodiment shown in Figures 1 
and 2, the sensor comprises two thin film members 
32 and 34, member 32 comprising sensor 22 and 34 
comprising sensor 24, each member comprising 
one-half of heater 26 and having a preferred 

60 dimension of 150 microns wide and 400 microns 
long. 

The embodiment of the system further describes 
an accurately defined air space 30 which contem- 
plates air space effectively surrounding elements 22, 
65 24, 26. The effectively surrounding air space is 
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achieved by fabricating the structure on silicon 
surface 36, thin film elements 22, 24 and 26 having a 
preferred thickness of approximately 0.08 to 0.12 
micron with lines on the order of 5 microns wide and 
spaces between lines on the order of 5 microns, the 
elements encapsulated in a thin silicon nitride film 
preferably having a total thickness of approximately 
0.8 microns or less, and by subsequently etching an 
accurately defined air space, of about 100 microns 
deep, into silicon body 20 beneath members 32 and 
34. 

Members 32 and 34 connected to top surface 36 
of semiconductor body 20 at one or more edges of 
depression of air space 30. As illustrated in Figure 3, 
members 32 and 34 may be bridged across 
depression 30; alternately, for example, members 32 
and 34 could be cantilevered over depression 30. 

Heat flows from the heater to the sensor by means 
of both solid and fluid couplings there between. Of 
note is the fact that silicon nitride (Si3N4) is a highly 
effective solid thermal insulator. Because the con- 
necting silicon nitride film within members 32 and 34 
is a good insulator, heat transmission through the 
solid does not dominate the propagation of heat 
from heater 26. This further enhances the relative 
amount of the heat conducted to sensing resistor 22 
and 24 from heater resistor 26 by flow through the 
surrounding fluid rather than through the supporting 
nitride film. Moreover, the supporting silicon nitride 
film has a low enough thermal conductivity that 
sensing resistor grids 22 and 24 can be located 
immediately adjacent or juxtaposed to heating 
resistor grid 26. Thus, sensing resistor grids 22 and 
24 are in effect suspended rigidly in the air space 
proximate heater resistor 26 and act as thermal 
probes to measure the temperature of the air near 
and in the plane of heater resistor grid 26. 

The operation of the system in sensing air flow is 
described in detail in the above-referenced U.S. 
patent 4,501,144. Typical circuit implementation is 
discussed briefly with reference to Figures 4 and 5 to 
add some insight. The heater control circuit illus- 
trated in Figure 4 uses a Wheatstone bridge 46 
which further typically includes heater resistor 26 
and a resistor 40 in its first leg and a resistor 42, heat 
sink resistor 38, and a resistor 44 in its second leg. 
An error integrator includes amplifiers 48 and 50 
keeps bridge 46 balanced by varying the potential 
across it and thus the power dissipated in heater 
resistors 26. 

The circuitry of Figure 5 monitors the resistance 
difference between downstream sensor 24 and 
upstream sensor 22. This circuitry includes a 
constant current source 52 comprising an amplifier 
72 and a differential amplifier 54 further including 
amplifiers 68 and 70. The constant current source 
drives a Wheatstone bridge comprising two high 
impedance resistors 56 and 58 in one leg and the 
two sensing resistors 22 and 24 with a nulling 
potentiometer 60 in the other leg. The gain of 
differential amplifier 54 is adjusted by potentiometer 
62. Output 64 provides an output voltage that is 
proportional to the resistance difference between 
the two sensing resistors 22 and 24. 

To get some concept of the small size of the 



mlcrobridge, the power required by heater resistor 
to heat such a device 200° C, for example, above 
ambient temperature Is less than 0.010 watt. The 
exceedingly small thermal mass of the heater and 

5 sensor element structures, their excellent coupling 
to the surrounding fluid because of a high surface/ 
volume ratio, and the thermal insulation provided by 
the thin silicon nitride connecting them to the 
supporting silicon body, and the surrounding air 

10 space, all contribute to produce a system well suited 
to fast and accurate sensing. Response time 
constants as short as 0.005 second have been 
measured. Consequently, sensor elements can 
respond very rapidly to proximate environmental 

15 changes. 

Now with reference to the implementation of the 
present invention, Figures 7a, 7b, and 7c, depict 
three slightly differing embodiments or configura- 
tions representative in terms of number and arrange- 

20 ment of the heaters and sensors which can be used 
in this invention, in Figure 7a, in contrast to Figure 1, 
all of the elements 122, 124 and 126 are used as 
heaters. Figure 7b is an embodiment which is similar 
to the embodiment of Figure 1 with thin film element 

25 126 acting as heater and elements 122 and 124 
acting as sensors. The embodiment of Figure 7c, 
represents the preferred arrangement in which the 
element 122 acts as heater and element 124 acts as 
sensor. The effective gap and thus the thermal 

30 isolation between heater and sensor is desirably 
wider in the embodiment of Figure 7c. 

The actual general geometric structure of the 
embodiments of Figures 1-3, and 7a-7c is more 
clearly illustrated in the scanning electron micro- 

35 graph (SEM) photo of Figure 8. The precision with 
which the cavity and bridge elements are defined 
and located in spaced relation, as Figure 8 depicts, is 
particularly noteworthy. The SEM represents a 
magnification such that the inidated length of 0.010" 

40 appears as shown. 

In the implementation of the invention disclosed 
herein particular attention is directed to (1) setting 
specific temperature markers in the sensor to 
determine the time periods needed for achieving the 

45 corresponding temperature changes, (2) using tem- 
perature sensors which are physically separated 
from the heater so that the direct influence of the 
heater and heat conducted to the sensor other than 
via the fluid of interest is reduced, and (3) using a 

50 pulse which reaches at least a momentary steady- 
state plateau to determine k, which then is used with 
the transient measure to determine c p . 

Figure 6 graphically depicts a square wave 
electrical energy pulse 130 to the heater as at 126 

55 which results in quasi square wave heat pulses 
released by the heater. These in turn, result in 
reactive curves as at 131, 132 and 133 at the sensor 
which vary as described below. The pulse applied to 
the heater, for example, may have a height of about 4 

60 volts with a pulse width of 1 00 ms. Since the heater is 
closely coupled through the fluid medium to the 
sensors, the family of curves 131, 132 and 133 
resembles the shape of the input pulse 130. They 
show the heat response in the sensors 122 and 124. 

65 Figure 12 is a photograph of one oscilloscope trace 
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showing temperature rise and fall versus time for dry 
air at atmospheric pressure. It uses a different scale 
for time than does Figure 6, but illustrates the curve 
form produced by the pulsed input. The curves 
generally include beginning and ending transient 
portions flanking a relatively steady-state central 
portion. The relatively quick response of the sensor 
allows a relatively long steady-state to exist even 
with a pulse of 100 ms. Of course, the curves are 
affected by factors such as pressure and tempera- 
ture as they influence the effective thermal conduc- 
tivity and specific heat of the particular fluid of 
interest. 

Heat flowing from the heater element or elements 
to the sensor element or elements is conducted 
both through the fluid and through the solid 
semiconductor element support substrate or the 
like. It is advantageous with respect to the measure- 
ment of k or c p of the fluid of interest that the amount 
of heat reaching the sensor through the solid 
connections be minimized so that substantially all 
the measured thermal effect is generated via the 
fluid of interest. 

With respect to the transfer of heat to the 
sensor(s) some background information regarding 
the propagation of heat or temperature waves is 
presented. The speed of propagation, v t of a one 
dimensional wave (if it features an exponential decay 
profile) is constant and given by the expression: 
v - Dr/a - (Dr/b) 05 , (1) 
where: 

a is an exponential decay constant 

b is the rise time constant at a fixed location and 

Dt is the thermal diffusivity. 

A complete list of nomenclature and subscripts 
with units appears in Table I, below. Dr is related to k 
and c p by the expression Dt *» k/c P (2) 

Dt, therefore, if known, may be a key to obtaining 
c P . The rise time constant, b, was measured to be 
about 4 msec. For typical gases, Dr ranges from 1.7 
cm 2 /s for He to .054 cm 2 /s for C3H8. Metals exhibit 
high values such as 1.7, 1.1 and .18 cm 2 /s 
respectively for Ag, Cu and Fe. Insulators, however, 
are even lower than the gases at .004 cm 2 /s for glass 
and .0068 cm 2 for Si3N4 which, as discussed above, 
is a good insulator. The propagation speed, v, in a 
typical gas sample then is about (1/0.004) 05 - 15 
cm/s. This compares with (0.0068/0.004)° 5 « 1.3 
cm/s for Si3N4, assuming that the same rise time 
constant of about 4 ms is applicable to both the one 
measured in the &3N4 and the actual one in the gas. 

The effect is that the influence of the temperature 
wave propagating from one thin film strip, that is, the 
heater, to a second thin film strip, the sensor, both 
being embedded in a membrane of &3N4, is faster 
for the gas than for the Si3N4. This also supports the 
choice of a material such as Si3N4, since it reduces 
the contribution of heat flow through the solid media. 
This is beneficial to the accuracy of the system. 

Typical microbridge embodiments are illustrated 
by Figures 7a - 7c. They will now be explained in 
greater detail. 



TABLE I 



NOMENCLATURE 



5 



20 



oymooi 




Units 


a 


Exponential Decay Constant 


cm 


ai-an 


Constant 




A 


Area of Heat Transfer to 


cm 2 




Microbridge or to Gas 




b 


Rise Time Constant at a 


°C/s 




Fixed Location 




c p 


Specific Heat 


cal/(cm 3 °C) 


Dt 


Thermal Diffusivity, 


cm 2 /s 




Dt « k/Cp 




k 


Thermal Conductivity 


cal/(sm to C) 


L 


Length of Thermal 


cm 




Conductance Path in Gas or 






Solid 




P 


Pressure of Gas 


psia 


Q 


Power of Heat Release Rate 


watts 


Ro 


Resistance at Room 


ohms 




Temperature 




t 


Time 


s 


T 


Absolute Temperature 


°C 


U 


Bridge Output or Amplified 


V 




Bridge Output 




V 


Volume of Gas or Solid 


cm 3 




(Microbridge) 




V 


Speed of Propagation 


cm/s 


x 


Temperature coefficient of 


°C-1 




resistance 





SUBSCRIPTS 

35 c Conduction 

S Microbridge or Solid 

g Gas 

o Room, Reference or Gas 
Temperature Without 

40 Microbridge Heating 

h Heater or Hot 

m Middle or Medium 

45 The configuration of Figure 7a involves using the 
same microresistance 122, 124, 126 for the heating 
pulse and the sensing task. In this embodiment of 
the resistive heater-sensor element may be one leg 
of a conventional resistive Wheatstone bridge in a 

50 control circuit. 

Figure 7b depicts an arrangement wherein the 
center microresistance structure 126 is used as a 
heater flanked by two symmetrically located outer 
sensing resistance elements 122 and 124. The 

55 elements 1 22 and 1 24 are separated from the heater 
126 by a narrow gap. 

Figure 7(c) shows an embodiment configuration 
in which the left element of the bridge 122 is used as 
the heating element and the right element 124 as the 

60 sensor. This embodiment takes advantage of a 
rather large central gap to achieve improved thermal 
isolation between the heater and the sensor. 

Figure 9 shows a modified control circuit which 
uses the center microresistance 126 as heater, while 

65 the sensing task is performed by the two resistors 
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122 and 124. The dual heater sensor configuration 
corresponds to Figure 7b and the circuit is rep- 
resentative of typical sensor/measurement circuit. 
Figure 9 includes a timer 140 providing square-wave 
electrical pulses to the heater 126. The heater 
couples the heat pulse to the sensors 122 and 124 in 
the bridge 142. The output of the bridge is 
connected through an amplifier 143 to a pair of 
comparators 144 and 145 which operate "start" and 
"stop" inputs to a counter 146 which counts 10 mHz 
clock pulses. The counter counts measure the time 
interval (t2-ti) between temperatures T2 & T1 
illustrated in Figure 6. 

Figure 9a is similar to Figure 9 f but more detailed. 
The bridge configuration is the heater -space-sen- 
sor configuration of Figure 7c. The sensor resist- 
ance arm of the microbrldge is set into a Wheatstone 
bridge 150 at 124. Another proximate resistive arm 
122 is fed a voltage pulse from pulse generator 161 
to provide a heat pulse into the microbridge element 
126. The Wheatstone bridge 150 also may contain a 
nulling balancing resistor 152 which can be used in 
the manner of potentiometer 60 in Figure 5 to initially 
zero the device. The microbridge resistor sensor 
124 in the Wheatstone bridge receives the heat 
pulse from heater element 122 principally by thermal 
conduction through the surrounding fluid. Some 
conduction, of course, does occur through the solid 
microbridge substrate and surroundings. 

The circuitry of Figure 9a is conventional and can 
readily be explained with reference to its functional 
operation with regard to processing the bridge 
output signal. The voltage output signals of the 
bridge 150 are amplified by differential amplifiers 153 
and 154 in a differential amplifier section. The 
imbalance signal is further amplified by a high gain 
amplifier at 155. The signal at 156 as is the case with 
the signal at 147 in Figure 9 is in the form of a DC 
voltage signal, U, the amplitude of which is solely 
related to the thermal conductivity of the fluid of 
interest as will be discussed above. 

The remainder of the circuitry of Figure 9a 
includes a DC level clamping amplifier 157 and 
isolation amplifier 158. The temperature level, time- 
related switching and counting circuitry includes 
comparators 159 and 160 together with Nand gates 
161 and 162 having outputs which are connected to 
the counter timing device (not shown) as in Figure 9. 
By measuring the time needed for the sensor 
temperature to rise or fall between two or more 
known temperature values or markers as rep- 
resented by sensor resistance or bridge voltage 
outputs a measure related to the specific heat per 
unit volume, c p of the fluid of Interest is obtained. 
The timing device may be a conventional 10 MHz 
pulse counter or the like. Again, this is illustrated 
schematically in Figure 6. 

The output signal from the Wheatstone bridge, U, 
represents the voltage imbalance caused by the 
temperature change in microbridge sensor or 
sensors induced by the corresponding heater pulse 
output. Because the magnitude of this imbalance is 
related directly to the amount of energy absorbed by 
the sensor or sensors, the amplitude of the signal is 
directly related to the thermal conductivity, k, of the 



conducting media in a manner next explained. 

Figure 6 shows that during much of the about 
100ms wide pulse period the temperature of the 
sensor reaches and maintains a constant value. 

5 During this time, the influence of the energy sink or 
source terms represented by specific heat are zero, 
which means that only thermal conductivity governs 
the value of the sensor temperature. 
Figure 12 is a plot of temperature rise inthe form 

10 of bridge output, U, (Figure 9 or 9a) using the 
sensing arrangement of Figure 7(b) versus time in 
milliseconds for various gases at atmospheric 
pressure. Curves for methane, dry air, ethane and a 
vacuum are presented. In this specific embodiment 

15 there was a heater resistance of 800 ohms, a pulse 
height of 2.5 volts, and a pulse width of 100 ms. 
Temperature markers t, and X2 are shown on the 
graph. These markers relate to those of Figure 13 
which shows a graphical presentation of heat up 

20 time versus pressure for several gases with a 
sensor-heater such as that shown in Figure 7b and 
using the T2-T1, marked in Figure 11. 

The literature value of the thermal conductivity of 
several gases has been plotted vs. the measured 

25 sensor temperature expressed directly in terms of 
the measured Wheatstone bridge imbalance poten- 
tial, U. This relationship has been derived empirically 
for a microbridge of the type depicted In Figure 7(c) 
and Is plotted in Figure 13, using the least squares 

30 method in a multiple regression analysis to achieve 
the best fit curve. The relation can be linearized over 
a modest span sufficient for the purpose of the 
invention. Other combination configurations of 
heater/sensor embodiments can likewise be cali- 

35 brated using known gases or gases of known k. 
Thus, using an off-the-shelf flow sensor of the type 
7(c) in the circuit 9(a), a 4.0V pulse of 100 ms 
duration was used. 
This yielded an approximate linear relationship 

40 between U and kg of the form 
kg - a 4 U + as (3) 
where 

a4 = -25.8807 and as - 181.778 for the above 
conditions. 

45 The above then achieves the calibration of the 
sensor for k g . The linear approximation holds over 
enough of a span to provide accurate measure- 
ments. Similar relations may be derived under other 
measurement conditions including additional press- 

50 ure correction terms. 

Further details related to determinng the coeffi- 
cients for the algorithms to compute c p are 
described next. This determination requires that the 
measuring system be calibrated first, which consists 

55 of determining the coefficients ai, a2, and a3, of the 
alogirthm to then computer Cp. 

Assuming a two-dimensional model for heat 
transfer in the microbridge, see Figures 7a-7c, the 
measured sensor temperature response may be 

60 described with reference to the following processes 
(at zero gas flow): 

1 ) Heat release by the heater element film. 

2) Temperature build up in the heater 
element material (FeNi or Pt) and surrounding 

65 support material (Insulator Si3N 4 ), i.e. within the 



6 



11 



EP 0348 245 A2 



12 



bridge material. 

3) Conduction towards the sensor via a) the 
bridge material, and b) the fluid phase sur- 
rounding the bridge. 

4) Temperature build up in the sensor 
material (as in heater material in item 2 above), 
and In the gas surrounding it by the heat 
arriving via the above processes. 

5) Achieving a steady-rate distribution of 
temperature. 

6) The revenue process to steps 1-5 during 
the start of the heater off-period. 

Further assuming, for the sake of simplicity, that 
the specific heats of the involved gaseous and solid 
materials do not depend on temperature, we can 
approximately describe the above processes by the 
following expressions (see Table I above for symbol 
explanation) using the same process numbering as 
above: 

1) q =, V 2 /(R 0 (1 + a (Th-To)) for small 
temperature rises. 

2) The heater temperature results from 
balancing the heat input and output rates: 
Th-To - Q/(k»A*/U + Mo/U) with Q in watts; 
the temperature Th is established in a time that 
is short compared to the time it takes to reach 
the sensor if the sensor is not identical to the 
heater, as in configurations 7(b) and 7(c). 

3) In a truly one-dimensional case most of 
500/o of the released power Q eventually arrives 
at the sensor, since it only has two ways to go 
( + x and -x directions). In a two- (or even 
three-) dimensional case a major part of Q gets 
dissipated in the y and z directions, so that only 
a fraction, Q c , is conducted to the sensor, with a 
corresponding drop of the original temperature, 
Th, down to an intermediate temperature T m . 
The sensor then experiences an energy rate 
arrival of 

Qc - (Tm-To) (ksAs/U + kgAfl/Lg) (4) 

4) The sensor temperature rise rate is 
governed by the specific heat of the gas 
surrounding the sensor and the closely coupled 
material of the sensor itself so that: 

Q c = (dT/dt) Cp»V 8 + (dT/dt)CpgV fl (5) 
The quantity measured and plotted in Figures 14, 
15 and 16, is the time (dt) needed to raise the sensor 
temperature by an increment (dT) which is chosen 
by the two or more sensor resistance value markers 
corresponding to Tn and T2. 

It is readily apparent from equation (5) that c w 
could be determined for an unknown gas if the 
various quantities entering in Eqs. (4) and (5) were 
either known or measurable. It has been found, 
however, that even if only dt, dT, T 0 , P and k g are 
conveniently measurable, the other quantities may 
be determined by calibration. This can be done 
according to an invention as follows: 

For calibration, gases of known composition 
(preferably but not necessarily pure) and therefore 
of known specific heat and thermal conductivity at 
the used pressure and temperature (both also 
measured), are brought in contact with the sensor. 
The effect of the pulsed heat releases is recorded in 
terms of the lapsed time, t2-ti, as has been 



described. After noting results for various gases, 
pressures, heater temperatures and/or heating/ 
cooling periods, with pulses of constant tempera- 
ture, voltage, current or power, the recorded time 
5 and condition data are entered into an array of data 
ports which can be used for automatic or compu- 
terized data processing or other number crunching 
techniques. 

The process can be illustrated with the help of 
10 equations (4) and (5), by way of example, without 
excluding other, similar approaches likely to occur 
to one skilled in numerical analysis. With this in mind, 
the following ports receive data or input for various 
gases, pressures (and temperatures): 

15 

Ports: Y X1 X2 

Inputs: CpgP/Po (t2-ti)k g ts-ti 

20 Known and available multiple linear regression 
analysis (MLR A, see Figure 10) program can 
determine the linear coefficients ai, a2, and a3 (e.g., 
by matrix inversion), which, together with the above 
input data, forms the calibrated expression derived 

25 from equations (4) and (5) to compute specific heat, 
c P :c pg P/Po = ai(t 2 -ti)k g + a 2 (t 2 -ti) -as (6) 

The determined (calibration)coefficients, of 
course, represent the lumped factors of several 
sensor properties or conditions from equations (6) 

30 and (7): 

ai = (Tm-To) (A fl /L g )/(V g dT), 

a 2 - (T m -To) (Ao/Ls)/(V g dT)k s , (7) 

a3 - CpsVs/Vg 

In order to minimize differences in T m at the 

35 sensor location, the most advantageous operation 
from among constant temperature, voltage, current 
or power is chosen. The above method is demon- 
strated on the basis of 1) constant voltage pulses, 
which result in quasi square wave heat pulses 

40 released by the heater, and 2) changes in gas type 
(CH4, C2H6, air and O2) and pressure; the chosen 
configuration was 7(b). 

Figure 14 shows the result of storing and plotting 
the dt = t2-ti and pressure data for each of the 

45 gases used, for which the c P and k values can be 
obtained from the open literature. This relation is 
linearized by applying the least squares method in a 
multiple linear regression analysis to achieve the 
best fit line. After entering these data into the above 

50 ports Y, X1 and X2, the regression analysis program 
performed. The obtained result was, for a configura- 
tion as in Figure 7(b): 

ai = -16509, a 2 - 3.5184 and a 3 = .005392 (7a) 
Proof that the above calibration coefficients are 

55 valid is provided by Figure 15, for example, in which 
these coefficients have been used to generate the 
shown lines for CH 4 , C2H6, air and O2. As shown, the 
lines indeed connect and agree with all experimental 
points. Additional lines have been plotted with the c P 

60 and k data of the literature for other gases as well. 
The final step in using this calibration method 
involves known means to store, write or burn in the 
obtained, tailored values of ai, a2 and a3 for the 
individual microbridge, which may be a Honeywell 

65 MICRO-SWITCH Model No. AWM-2100V, into the 
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memory linked to It. The microsensor Is then ready 
for use to measure the specific heat of unknown 
gases, provided that P and k be known at the time of 
measurement. 

Figure 10 depicts a schematic block diagram of a 
device for measuring Cp and k. The system Includes 
the signal processing circuitry indicated by 170, a 
multiple linear regression analysis (MLRA) unit 171 
for deriving the known equation constants for the 
particular microbridge configuration and circuitry 
used, i.e., ai - an, a data bank 72 for storing 
calibration c P and k data and an output interface unit 
173. 

With respect to the embodiment of Figure 10, 
prior to use, field recalibration may be accomplished 
simply by entering the P, c P and k values of the test 
gas into the data bank. If P cannot be measured 
independently of the sensor already In the subject 
system its errors can be incorporated as a correc- 
tion in the c p and k recalibration. The measured 
values of U and dt are then used as in the 
measurement mode to determine sensor values of k 
and Cp. If they disagree from the entered values the 
constants az and as may be modified to fit the 
entered or book values. 

This approach may be a practical one for field use, 
but it should be checked by using a second test gas. 
If that agrees, the recalibration may be completed. If 
not, a complete calibration of all ai-as coefficients 
should be made. 

It should be mentioned that in all of the above 
discussion the influence of temperature was not 
mentioned for the sake of simplicity. It is well known, 
however, that temperature does influence both c p 
and k but can be addressed, if necessary, in one of 
the following ways: 

1) Controlled, (expensive and energy con- 
suming) or 

2) Compensated by special temperature- 
sensitive elements in the analog part of the 
circuit, or 

3) Entered into the sensor algorithm as an 
additional parameter, which is sensed, e.g., by 
monitoring one of the many available tempera- 
ture dependent resistors on the sensor. This is 
the preferred approach for sensing systems 
requiring maximum accuracy. 

With respect to use of the instrument of Figure 10, 
the U and dt = t2-ti (and P) signals obtained for an 
unknown gas are processed as follows in this mode: 

1) Computation of k from expression (3) 
using the coefficients a 4 and as which have 
been stored in (or burned into) the sensor's 
memory after calibration, and 

2) Computation of c p from expression (6). It 
should also be noted that a pressure signal is 
also needed as a basic ingredient since c p is 
used here in relation to a volume of gas as 
opposed to k which is largely pressure inde- 
pendent if the sensor is used at or above 
atmospheric pressure, at which the gas mean 
free path is small compared to the characteris- 
tic dimensions of the involved sensor. 

The graphical presentation of Figure 16 depicts 
heating time in milliseconds versus pressure and 



gas type and specifically showing curves for meth- 
ane, ethane, air and oxygen. The sensing configura- 
tion of Fibure 7(c) was used. In this example, the 
pulse height was 1 .75 volts with a pulse width of 100 

5 ms. and the heater and sensor resistance each 
being about 2000 ohms. Figure 17 depicts a cooling 
curve for the same configuration as Figure 16. 
Conditions were the same except that the pulse 
height was 4.0 volts. 

10 Of course, the output of the device can be in any 
desired form including analog or digital signals, 
printed records, etc., after the value is obtained. 



15 Claims 

1. A method for determining the thermal 
conductivity, k, and specific heat, c p , of a fluid of 
interest using proximately positioned hearing 

20 (126,152) and sensing means (122,124) coupled 

by said fluid of interest, said sensing means 
(122,124) having a temperature sensitive out- 
put, characterised by the steps of: 
providing an energy input to the heater means 

25 (126,152) of a level a*nd of a duration such that 

both transient temperature change and a 
substantially steady-state temperature occur in 
the sensing means; 

determining k of the fluid of interest based upon 
30 a known relation between the sensor output 

and k at steady-state sensor temperature; and 
determining c P of the fluid of interest based on a 
known relation among k, the rate of change of 
sensor output during a transient temperature 
35 change in the sensor and c p . 

2. The method according to Claim 1 charac- 
terised in that the heater and sensor means are 
electric resistance elements and the input to 
the heater is in the form of an electrical pulse. 

40 3. A method for determining the thermal 

conductivity, k, and specific heat, c P , of a gas of 
interest using proximately positioned micro- 
bridge electrical resistance heating (126,152) 
and thermal sensing means coupled by said 

45 fluid of interest, said sensing means (122,124) 

having a temperature sensitive electrical output 
signal characterised by the steps of: 
providing an electrical energy input pulse to the 
heater means (126,152) of a level such that the 

50 sensing means (122,124) experience a transient 

temperature change and of a duration such that 
a substantially steady-state temperature is 
achieved in the sensing means; 
determining k of the gas of interest based upon 

55 a known relation between the electric sensor 

output signal and k at steady-state sensor 
temperature; and 

determining c p of the gas of interest based on a 
known relation among k, the rate of change of 

60 sensor output during a transient temperature 

change in the sensor and c p . 

4. A method for determining the thermal 
conductivity, k, of a fluid of interest using 
proximately positioned microbridge electrical 

65 resistance heating (126,152) and thermal sens- 
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ing means (122,124) coupled by said fluid of 
interest, said sensing means having a tempera- 
ture sensitive output characterised by the steps 
of: 

providing an electrical energy input pulse to the 5 
heater means of a known level such that the 
thermal sensing means experiences a transient 
temperature change and of a known duration 
such that a substantially steady-state tempera- 
ture is achieved in the sensing means and ; 10 
determining k of the fluid of interest based upon 
a relation between the sensor output and k at 
steady-state sensor temperature substantially 
approximated by 

k « a4U + a 5 15 
where 

U is the sensor output and 
a* and as are constants. 

5. A method for determining the thermal 
conductivity, k, and specific heat, c P . of a fluid of 20 
interest using proximately positioned micro- 
bridge electrical resistance heating (126,152) 

and thermal sensing means (122,124) coupled 
by said fluid of interest, said sensing means 
having a temperature sensitive output charac- 25 
terised by the steps of: 

providing an energy input pulse to the heater 
means of a known level such that the sensing 
means experiences a transient temperature 
change and of a known duration such that a 30 
substantially steady-state temperature is 
achieved in the sensing means; 
determining k of the fluid of interest based upon 
the relation 

k = a 4 U + as & 
where 

U is the sensor output at steady state and 

a4 and as are constants; and 

determining c P of the fluid of interest based on 

the relation among k, the rate of change of 40 

sensor output during a transient temperature 

change in the sensor and c p 

CpP/Po = ai(t 2 -ti)k + a 2 (t 2 -ti)-a 3 

where 

ai , B2 and a3 are constants 45 

P = pressure (psia) 

P 0 = reference pressure (psia) 

(t2-ti) = measured time span between known 

temperatures. 

6. The method according to any one of the 50 
preceding Claims characterised by the further 

step of providing an output signal indicative of 
the value of k and/or c p . 

7. The method according to Claims 3 or 5 
characterised in that c p is determined with 55 
respect to an upward transient temperature 
change in the sensor. 

8. The method according to Claim 3 to 7, 
characterised in that c p is determined with 
respect to a downward transient temperature 60 
change in the sensor. 

9. Apparatus for determining thermal con- 
ductivity, k, and specific heat. c p , of a fluid of 
interest characterised by; 

heat means (126,152); 65 



thermal sensing means (122,124) in proximate 
position to said heater means and in communi- 
cation therewith through the fluid of interest, 
said sensing means characterised by a tem- 
perature dependent output; 
energizing means (140,151) connected to said 
heater means (126,152) for energizing said 
heater on a time-variable basis in a manner to 
induce both transient and substantially steady- 
state temperature conditions in said thermal 
sensing means; 

first output means (143,155) for providing first 
output signal (156) indicative of the tempera- 
ture of said thermal sensing means; 
means (144-146,158-162) for determining the 
rate of change of temperature of said tempera- 
ture sensing means; 

means (170-173) for determining k of the fluid of 
interest based upon a known relation between 
the first output and k at steady-state sensor 
temperature; and 

means (170-173) for determining c p of the fluid 
of interest based on a known relation among k, 
the rate of change of the first output during a 
transient temperature condition and c p . 

10. Apparatus for determining the thermal 
conductivity, k, and the specific heat, c p , of a 
fluid of interest, characterised by: 
amicrobridge system (142,150) including resis- 
tive heater portion (126,152) and resistive 
sensor portion (122,124) in juxtaposed spaced 
relation, said heater and sensor portions each 
having terminals, said system further being 
positioned in direct communication with the 
fluid of interest said resistive heater thereby 
being thermally coupled to said sensor via said 
fluid of interest; 

electrical pulse producing means (140,151) 
connected in energizing relation to said heater 
terminals for providing an energy input to the 
heater of a level and duration such that both a 
transient and substantially steady-state tem- 
perature conditions occur in the sensing 
means; 

first output means (143,155) for providing an 
electrical potential output signal (156) indicative 
of the temperature of said sensing means; 
means (144-146, 158-162) for determining the 
rate of change of temperature of said sensing 
means; 

means (176-179) for determining k of the fluid of 
interest based upon a know relation between 
the sensor output and k at steady-state sensor 
temperature; and 

means (170-173) for determining c p of the fluid 
of interest based on a known relation among k, 
the rate of change of sensor output during a 
transient temperature change in the sensor and 
c p . 

11. The apparatus according to Claim 9 or 10 
characterised by second output means for 
providing an output indicative of k and/or c P of 
the fluid of interest 

12. The apparatus according to any one of 
Claims 9 to 11 characterised in that said 
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THERMAL DIFFUSIVITY SENSOR OR 
SPECIFIC HEAT SENSOR 
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P, PRESSURE IN PSIA 

Fig. 17 

COOLING TIME VS. PRESSURE AND GAS TYPE, 
USING AN OFF-THE-SHELF FLOW SENSOR 
OF CONFIGURATION 7b 
SENS-HTR-GAP-HTR-SENS 
PULSE HIEGHT: 2.5V 
WIDTH: 100ms 
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